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Abstract The properties of the first galaxies are shaped in large part by the first generations of stars, which emit 
high energy radiation and unleash both large amounts of mechanical energy and the first heavy elements when they 
explode as supernovae. We survey the theory of the formation of the first galaxies in this context, focusing on the 
results of cosmological simulations to illustrate a number of the key processes that define their properties. We first 
discuss the evolution of the primordial gas as it is incorporated into the earliest galaxies under the influence of the 
high energy radiation emitted by the earhest stars; we then turn to consider how the injection of heavy elements by the 
first supernovae transforms the evolution of the primordial gas and alters the character of the first galaxies. Finally, we 
discuss the prospects for the detection of the first galaxies by future observational missions, in particular focusing on 
the possibiUty that primordial star-forming galaxies may be uncovered. 



1 Introduction: Defining Characteristics of the First Galaxies 

While the first stars are for the most part well-defined objects, the definition of the first galaxies is somewhat more 
ambiguous (see e.g. Bromm & Yoshida 2011). Here we shall adopt the common view that a galaxy must be able to 
host ongoing star formation, even in the face of the radiative and mechanical feedback that accompanies the formation 
and evolution of stars. By this definition, the formation sites of the first stars, dark matter minihalos with masses 10^ 
- 10*' M,.„ are unlikely candidates for the first galaxies, as the high energy radiation emitted by young stars and the 
supernovae that mark their end of life can rarify and expel any dense gas from which stars may form at a later time. As 
shown in Figure 1, it is only somewhat larger halos, with masses 10^ - 10^ M©, which have deep enough gravitational 
potential wells and enough mass to prevent the expulsion of gas after an episode of star formation (e.g. Kitayama & 
Yoshida 2005; Read et al. 2006; Whalen et al. 2008). 

As can be inferred from this Figure, one of the distiguishing characteristics of halos massive enough to host ongoing 
star formation, and so to host the first galaxies, is the characteristic temperature Tyir that gas reaches during their 
virialization. This, referred to as the virial temperature of the halo, can be derived by assuming that the absolute 
magnitude of the gravitational potential energy of the halo is twice its kinetic energy, which yields 

where Mh is the mass of the halo, z is the redshift at which it collapses, and /i is the mean molecular weight of the 
gas in the halo, here normalized to a value appropriate for neutral primordial gas. The Hubble constant = 100 h 
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Fig. 1 The fate of host halos of the first stars, 
as a function of the mass Mj, of the halo and 
the energy £sN with which the stars explode as 
supemovae. Less massive halos have their gas 
blown away when even relatively weak (e.g. 
£sN ~ 10^" 2rg) Pop III supemovae explode 
within them, both along with {circles) and 
in the absence of (triangles) the additional 
radiative feedback from the progenitor stars. 
Halos which retain their gas are shown by 
crosses. From equation (1) we can see that for 
the suite of halos shown here at z ~ 20, those 
with Tyir > 10* K are able to retain their gas, 
in contrast to the less massive minihalos. From 
Kitayama & Yoshida (2005). 
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km ^ Mpc^ ' also appears here through h. Q From Figure 1 we see that the mass of halos which are large enough to 
host ongoing star formation, at z ^ 20, is ^ 10^ M©; this corresponds to a virial temperature of ^ 10'* K. One 
of the reasons for this is that 10'* K is roughly the temperature to which photoionization by stars heats the gas (see 
e.g. Osterbrock & Ferland 2006); thus, gas that is photoheated by stars remains bound within a halo with such a virial 
temperature. In turn, the presence of this gas when stars explode as supernova leads to the rapid loss of the mechanical 
energy in the explosion to radiation, thereby limiting the amount of gas blown out of the halo, in contrast to the case of 
the first supemovae in less massive minihalos (see Section 3.1). Also, due to the efficient cooling of atomic hydrogen 
at this temperature, gas can collapse into halos with Tyir > 10"* K regardless of its molecular content, in contrast to the 
minihalos that host the first stars, into which primordial gas only collapses if it is cooled by H2 molecules (e.g. Oh & 
Haiman 2002); this implies that star formation can take place even under the influence of the molecule-dissociating 
radiation emitted by the first stars (see Section 2.2). 

Figure 2 shows the properties of an atomic cooling halcQ , in which a first galaxy would form, at z ^ 10 in a cosmo- 
logical simulation (see Greif et al. 2008). As shown here, much of the primordial gas that falls from the intergalactic 
medium (IGM) into the potential well of the halo is shock-heated to Tyir ^ 10'* K at a physical distance of ~ 1 kpc 
from the center of the halo. This distance corresponds to the virial radius rvh- of the halo, defined in general terms 
as the radius within which the average matter density is equal to the value at which virial equilibrium is established, 
which is ~ IStt^ times the mean matter density of the universe at the redshift z at which the halo forms (e.g. Barkana 
& Loeb 2001). For the standard ACDM cosmological model, this is given in physical units as 



where we have normalized to values of halo mass and redshift that are typical for atomic cooling halos hosting the first 
galaxies. Near the virial radius a large fraction of the gas is hot (> 500 K) and rotating about the center of the halo at 
nearly the circular velocity Vdrc of the halo (Greif et al. 2008), defined as the velocity with which a body must move 
in order to be centripetally supported against gravity at the virial radius: 

^''^^^- 20r Vr^Vkms-'. (3) 



Mh Y fl 
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' Note that this formula is derived assuming a standard CDM cosmological model in which /; ~ 0.7 (see e.g. Barkana & Loeb 2001); as 
such, this formula is valid at the high redshifts (i.e. z » 1) at which the tirst galaxies form, but must be modified at lower redshifts in 
order to account for a cosmological constant A. 

^ Because the primordial gas can cool via emission from atomic hydi'ogen and collapse into halos with Tyii- ~ 10* K, such halos are 
commonly referred to as 'atomic cooling' halos. 
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Fig. 2 The properties of the primordial gas collapsing into an atomic cooling dark matter halo at z —10. Shown are the hydrogen number 
density (left panel) and temperature (right panel), the dashed lines denoting the virial radius ryj,- at a distance of ~ 1 kpc. Note that most of 
the gas is accreted directly from the IGM and shock-heated to the virial temperature of Tvii- ~ K, although cold accretion also becomes 
important as soon as gas cools in filaments and flows towards the centre of the galaxy, such as through the streams coming from the left- 
and right-hand sides of the panels. In contrast to the minihalos in which the first stars form, a halo with a virial temperature T^,n- > 10"* K is 
massive enough and has a deep enough gravitational potential well to retain its gas even when stars formed within it explode as supemovae 
(see Figure 1). Hence, such halos are strong candidates for the formation sites of the first galaxies. From Greif et al. (2008). 



However, there is also a substantial portion of the infalling gas that falls to the center of the halo in cool, dense filaments 
and is not shock-heated to the virial temperature. These dense filaments feed cold gas into the central ~' 100 pc of the 
halo, contributing to the majority of the gas the temperature of which is < 500 K and which may collapse to form stars 
(Greif et al. 2008). 

While the atomic cooling halo shown in Figure 2 is a prime example of the type of halo in which the first galaxies 
likely formed, there are numerous physical effects that were not included in the cosmological simulation from which 
this halo was drawn, most notably the feedback effects of Population (Pop) III stars (see e.g. Wise & Abel 2008; 
Johnson et al. 2008; Greif et al. 2010; Whalen et al. 2010). The high energy radiation emitted by the first stars both 
ionizes the primordial gas and dissociates molecules, which are critical cooling agents. Also, many of the first stars 
explode as violent supemovae, which inject large amounts of mechanical energy into their host minihalos and the 
IGM, as well as dispersing the first heavy elements, thereby altering forever the properties of the gas from which the 
first galaxies form. 

In this Chapter, we shall focus on how this feedback from the first generations of stars impacts the formation and 
evolution of the first galaxies. In Section 2, we briefly discuss how the cooling properties of the primordial gas, which 
shape the nature of Pop III star formation, are affected by the radiation emitted from the first stars and accreting black 
holes. In Section 3, we then turn to discuss how the first supemovae enrich the primordial gas with heavy elements, 
and how this process leads to the epoch of metal-enriched Pop II star formation. In Section 4, we briefly discuss the 
prospects for observing the first galaxies, and for finding Pop III star formation therein, using facilities such as the 
James Webb Space Telescope (JWST). Finally, in Section 5, we close with a summary of the results presented in this 
Chapter and give our concluding remarks. 
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2 Evolution of the Primordial Gas in the Formation of the First Galaxies 



Being composed solely of the hydrogen, helium, and trace amounts of lithium and beryllium synthesized in the Big 
Bang, the primordial gas contains a limited number of coolants, chief among these H2 at temperatures < 10"*^ K. 
Because of the ineffecient cooling of the gas relative to the metal-enrichec0 gas from which stars form today, it is 
Ukely that the Pop III initial mass function (IMF) is top-heavy compared to that of the stars observed in our Milky 
Way. A simple explanation for this is based on the mass scale at which the fragmentation of the primordial gas takes 
place. Known as the Jeans mass Mj, this is essentially the mass at which density enhancements grow via gravity more 
quickly than they can be erased due to pressure gradients. To estimate Mj for a gas with a number density n and a 
temperature T, related to the sound speed Cs by 2>kzTI2 = ^ninclll, we first estimate the timescale at which density 
enhancements grow as the free-fall time ~ (Gp)"5= (GlXmnn)^"!, where G is Newton's constant. Then, estimating 
the timescale in which density enhancements are erased as the sound-crossing time fsc — L/cs, we equate these two 
timescales to estimate the characteristic size Lj and mass of a gas cloud which is just massive enough to collapse under 
its own gravity. We thus arrive at an expression for the Jeans mass Mj, given by 



where we have assumed = 1.2, appropriate for neutral primordial gas, and have again normalized to quantities typical 
for primordial star-forming clouds. As we shall discuss in Section 3, the primordial gas is in general unable to cool 
as efficiently as metal-enriched gas, which leads in general to higher temperatures at fragmentation and so to a larger 
characteristic mass of the gravitationally unstable gas clouds from which stars form (e.g. Bromm & Larson 2004). 

While the Jeans mass is an estimate of the mass of a collapsing gas cloud, the amount of gas that is finally incor- 
porated into a star is also dictated by the rate at which gas accretes onto it, starting from the formation of a protostar. 
Thus, another reason that primordial stars are likely to be more massive than stars forming from metal-enriched gas is 
that higher gas temperatures also translate into higher accretion rates, as can be seen by estimating the accretion rate 
Mace as a function of the temperature of the gas (see e.g. Stahler et al. 1980). Assuming that, through the action of 
gravity, the protostar grows by accreting from a gas cloud of mass ~ Mj, the accretion rate can be estimated as 



where we have again assumed = 1.2 and normalized to the characteristic temperature of the gas from which Pop III 
stars form in minihalos (see e.g. Glover 2005). Therefore, it is a combination of both the relatively large reservoir of 
gas available in gravitationally unstable gas clouds and the relatively high accretion rates onto primordial protostars 
(e.g. Omukai & PaUa 2003; Tan & McKee 2004; Yoshida et al. 2008) which suggests that Pop III stars are more 
massive than metal-enriched Pop II or Pop I stars. As both Mj and Mace depend strongly on the temperature of the 
gas, one of the central questions with regard to star formation in the first galaxies is the degree to which the gas is 
able to cool. In the next Section, we discuss the cooling of primordial gas in the first galaxies, focusing on how it is 
different from the case of cooling in the minihalos in which the first stars form. Later, in Section 3.3, we discuss how 
the cooling properties of the primordial gas change when it is mixed with heavy elements and collapses to form the 
first Pop II stars. 



2.1 Cooling of the Primordial Gas 

As shown in Figures 2 and 3, the primordial gas collapsing into atomic cooling halos is typically shock-heated to 
the virial temperature of > 10'* K. In contrast to the case of Pop III star formation in minihalos, the gas at these 
temperatures is partially ionized, and this can have important consequences for the evolution of the gas as it collapses 



We use the common term 'metals' to refer to elements heavier than helium which are produced in stars and supemovae. 
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Fig. 3 The properties of the primordial gas inside a first star-forming minihalo at z = 23 {left panel) and a first galaxy-hosting halo at z = 
10 (right panel), as found in a cosmological hydrodynamics simulation. The temperature, electron fraction, HD fraction, and H2 fraction 
are shown as functions of number density, clockwise from top left to bottom left. Left panel: In the minihalo case, adiabatic collapse drives 
the temperature to T > 10^ K and the density to n > 1 cm^, where molecule formation sets in and allows the gas to cool to ~ 200 K. At this 
point, the central clump becomes Jeans-unstable and presumably collapses to form at least one massive Pop III star. Right panel: In the first 
galaxy, formed in an atomic cooling halo, a second cooling channel emerges due to an elevated electron fraction at the virial shock, where 
the temperature rises to ~ 10* K; this, in turn, enhances molecule formation and allows the gas to cool to the much lower temperatures, as 
shown by the evolutionary track of a representative parcel of gas (green lines). Due to the correspondingly lower Jeans mass (equation 4) 
and accretion rate (equation 5), less massive Pop III stars (Pop III. 2; see Section 2.1) are expected to form in such a first galaxy, perhaps 
with a characterisitic mass of the order of ~ 10 Mp,. From Greif et al. (2008). 



to form stars in the first galaxies. To see why, we note that the primary reaction sequence leading to the formation of 
H2 molecules is (e.g. GalU & Palla 1998; Glover 2005) 

e-+H^H- + 7 (6) 

H-+H^H2 + e-, (7) 

where 7 denotes the emission of a photon. Whereas the primordial gas which collapses into minihalos to form the 
first stars has a free electron fraction < 10^"*, the collisional ionization of the primordial gas collapsing into atomic 
cooling halos can lead to an enhancement of the free electron fraction by a factor of more than an order of magnitude, 
as shown in Figure 3. In turn, this leads to high rates of H2 formation in atomic cooling halos, principally via the above 
reactions for which free electrons act as catalysts (e.g. Shapiro & Kang 1987). The net result is a generally higher H2 
fraction in the high density, central regions of atomic cooling halos than in minihalos, as is also shown in Figure 3, 
and hence also to higher cooling rates due to molecular emission. Therefore, somewhat counter-intuitively, because of 
the higher virial temperatures of the atomic cooling halos in which the first galaxies form, the dense gas in the centers 
of these halos can cool more effectively than in the minihalos in which the first Pop III stars form. 

In fact, the ionization of the primordial gas in atomic cooling halos results in the formation of another molecule 
which can be even more effective at cooling the gas than H2: deutrerium hydride (HD). With the high H2 fraction that 
develops in partially ionized gas, HD forms rapidly via the following reaction; 
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D+ + H2^HD + H+. (8) 

While deuterium is less abundant in the primordial gas than hydrogen by a factor of the order of 10^^, the HD molecule 
is able to cool to temperatures considerably lower than H2 (e.g. Flower et al. 2000). Firstly, this owes to the fact that HD 
has a permanent dipole moment, allowing dipole rotational transitions, which spontaneously occur much more often 
than the quadrupole rotational transitions in H2. Also, the dipole moment of HD allows transitions between rotational 
states of Ay = ±1, which are of lower energy than the AJ = ±2 quadrupole transitions of H2. Thus, collisions with 
other particles, such as neutral hydrogen, can excite the HD molecule from the ground to the first excited rotational 
state (7=1), from which it decays back to the ground state by a dipole transition. The photon that is emitted in the 
process carries away energy and thus cools the gas. Because HD can be excited to the 7=1 state by relatively low 
energy collisions, and because its subsequent radiative decay occurs quickly compared to that of H2, the cooUng rate 
per molecule is higher for HD than H2 at temperatures < 100 K, as shown in Figure 4. 

Whereas in cosmological minihalos H2 cooling alone can cool the gas to ~ 200 K, as shown in the left panel of 
Figure 3, HD cooling can be so effective as to allow the primordial gas to cool to the lowest temperature that can be 
achieved via radiative cooling, that of the cosmic microwave background (CMB), Tcmb = 2.7(l+z) (e.g. Larson 2005; 
Johnson & Bromm 2006; Schneider & Omukai 2010). It is useful at this point to derive this fundamental result, as we 
will draw on the formalism introduced here later as well, in discussing the impact of the first heavy elements on the 
cooling of the primordial gas (see Section 3.3). 

To begin, note that the frequency Vio of emitted radiation for the rotational transition 7 = 1 of HD can be 
expressed as 

^^130K, (9) 

where is the Boltzmann constant and h is the Planck constant. For clarity, here we shall consider the simple case 
in which only this transition and its reverse occur Next, consider a finite parcel of primordial gas with a temperature 
Tgas- For simplicity, we shall assume that the density of the gas is sufficiently high to establish local thermodynamic 
equilibrium (LTE) level populations according to the Boltzmann distributioiQ 

ni gi - 



e *B?k- ^ (10) 
"0 go 

where «i is the number density of HD molecules in the /th excited rotational state and gi is the statistical weight of 
that state; specifically, here we have ^1 = 3^o- Furthermore, as we are considering only transitions between the ground 
state and the first excited state, we shall take it that no other rotational levels are occupied. Equivalently, we take it here 
that Tgas < 130 K, as otherwise collisions with other particles would be sufficiently energetic to excite the molecule to 
higher levels. Finally, we make the assumption that Tgas > Tcmb- Thus, if we denote the specific intensity of the CMB, 
which is an almost perfect blackbody, at the frequency Vio as /vm, then it follows that 

2/!vfo/c2 2/ivfo/c2 

Now, with the Einstein coefficients for spontaneous and stimulated emission from 7=1^0 denoted by A 10 and 
Bio, respectively, and that for absorption of a photon effecting the transition 7 = — > 1 by Bqi, we have the standard 
relations Bio^i = Boi^O ™d 

2/1 V?n All) 

= IT- ■ (12) 
Bio 

Along with equations (10) and (11), these imply that 



"* Due to infrequent particle collisions at low densities, the rate of radiative deexcitations can exceed that of coUisional deexcitations, leading 
to non-LTE level populations (see Section 3.3). 
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Fig. 4 Molecular cooling rates for H2 (solid 
line), HD (dot-dashed line), and (dashed 
lines). Owing primarily to the permanent 
dipole moment of HD, at low temperatures the 
cooling rate per HD molecule is significantly 
higher than that per H2 molecule. Thus, pri- 
mordial gas enriched in HD is able to cool to 
much lower temperatures than gas containing 
solely H2; indeed, as is shown in Section 2.1, 
primordial gas sufficiently enriched in HD 
can cool to the lowest temperature attainable 
by radiative cooling, that of the CMB. While 
the cooling rates per H2^ molecule (shown 
here for coUisional excitation by free elec- 
trons and hydrogen atoms) can be very high, 
the low abundance of this molecule limits its 
importance for the thermal evolution of the 
primordial gas. Finally, we note that at temper- 
atures ~ Iff* K, the primordial gas is cooled 
primarily by recombination and resonance 
lines of atomic hydrogen (not shown). From 
Yoshida et al. (2007a). 
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"O^Ol^Vio < «lAlO + «lBl0^v,o ■ (13) 

Thus, the gas is cooled, as more energy is emitted into the CMB radiation field than is absorbed from it. The rate at 
which the temperature drops can be found by first expressing the energy density of the gas as 



;nkBTg,, , (14) 



_ 3 

Mgas — 2' 

where « is the total number density of the gas particles, including all species. With this, equation (13) implies that, 
with no change in the density of the gas. 



/iVio[«oBoi/vio -"1^10 -niBio^vi,,] = ^n^B— ^ ■ (15) 



3 , dTsas 

Next, we take it that the ratio of the number density of HD molecules hhd to the total number density of particles « in 
the gas is given by the constant factor 

_ «HD «o + ni «o ,,^s 

ahd = — — — ■ (16) 

n n n 

Then using equations (10), (11), and (12) in equation (15), and neglecting stimulated emission for simplicity, the 
thermal evolution of the gas is approximately described by 

dt 

It is clear from this result that if Tcmb < T^as < 130 K, with the gas cooling only by radiative decay of the excited ro- 
tational state /=1 to /=0, the temperature of the gas will asymptotically approach Tcmb- Thus, equation (17) describes 
the fact that the CMB temperature is indeed a lower limit on the temperature to which a gas can cool via line emission 
only. Using the previous equation, we can estimate the timescale for reaching the CMB temperature floor as 

^ ' hvio 



e *b'cmb — e *B?gas _ (17) 



1 / ^bTcmb V^^^ / 



fcMB - ^, ^ — — ; exp 



(AioXhd)"' . (18) 
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Fig. 5 Expansion of the H 11 region around a 
Pop III.l star formed in a minilialo: tlie large- 
scale density distribution at j = 26 (top left), 
and the projected gas density at 1.1 (top right) 
and 3 Myr (bottom right) after the central 
star turns on. The bottom left panel shows 
the extent of the H II region (green) and that 
of the He III region (hhie) at 3 Myr When 
the central star turns off after this time, the 
strongly ionized primordial gas begins to cool 
and recombine, with H2 and HD molecules 
forming in abundance. This, in turn, enhances 
the cooling properties of the gas, which may 
later collapse and form Pop III. 2 stars in a 
more massive halo hosting a first galaxy. From 
Yoshida et al. (2007a). 
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Finally, we may use this timescale to define a critical HD abundance above which the gas may cool to the CMB, 
by demanding that the gas is able to cool faster than it is heated by compression during its collapse, which takes 
place roughly on the free-fall timescale. We thus require that fcMB ~ tff, where the free-fall time is calculated at the 
characteristic density n ^ 10^ cm^^ at which the primordial gas is found to fragment in cosmological simulations (e.g. 



Bromm & Larson 2004). With Ajo ~ 5 x lO^^ s" 
the critical HD abundance to be approximately 



for this transition (e.g. Nakamura & Umemura 2002), we thus find 
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-6 



(19) 



If the abundance of HD is lower than XHD,crit> the gas will not have time to cool to Tgas ~ Tcmb during its collapse. As 
shown in Figure 3, for the case of the primordial gas cooling in the first galaxies a large fraction of the gas at densities 
n > 10 cm^^ has an HD abundance greater than XnD.crit, whereas in the case of cooling in the minihalos hosting the 
first stars the HD fraction is in general much lower 

A high abundance of HD can, in general, be formed whenever the primordial gas becomes ionized. This occurs 
through collisional ionization, as in the case of shock heating to temperatures above ^ lO"* K in the virialization 
of atomic cooling halos, but also occurs when the first stars formed in minihalos emit high energy radiation which 
photoionizes the gas. As shown in Figure 5, a massive Pop 111 star emits enough ionizing radiation to destroy almost 
all of the neutral hydrogen within a distance of a few physical kiloparsec0 (e.g. Alvarez et al. 2006; Abel et al. 2007), 
via the reaction 



H + 7^H++e-. (20) 

Here the products are ionized hydrogen and a free electron, which has a kinetic energy equal to the energy of the 
ionizing photon minus the ionization potential of hydrogen, 13.6 eV. This free electron is ejected from the atom and 
shares it kinetic energy with other particles via collisions, thereby heating the gas to higher temperatures. Typically, an 
equilibrium temperature of ~ 10"* K is established in H II regions, largely set by a balance between the rate at which 
the gas is photoheated via the above reaction and the rate at which it is cooled by the radiative recombination and 
resonance emission of hydrogen (e.g. Osterbrock & Ferland 2006). While in the H II regions around active stars the 
temperature is thus too high for molecules to form in large abundances due to collisional dissociation, once the central 

' As H+ is also refeired to as H II, such photoionized regions fonued around stars are called H II regions. Likewise, the radiation froin 
massive stars, and especially massive Pop III stars, can doubly ionize helium within the so-called He III region (see Section 4). 
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star dies the hot ionized gas begins to cool and recombine. Under these conditions, molecules form rapidly and a high 
abundance of HD can be achieved (e.g. Nagakura & Omukai 2005; Johnson & Bromm 2006; Yoshida et al. 2007b; 
McGreer & Bryan 2008). 

Therefore, overall, primordial gas that has either been photoionized by a Pop III star in a minihalo or which has 
been partially ionized during the virialization of an atomic cooling halo, may in principle collapse and cool all the way 
to the temperature floor set by the CMB. This is distinct from the case of the first Pop III star formation in minihalos 
(e.g. Glover 2005), and this distinction motivates the following terminology (e.g. McKee & Tan 2008; Greif et al. 
2008; Bromm etal.2009il: 

Pop III. 1 The first generation of primordial stars formed in minihalos and not significantly affected by previous star 
formation. 

Pop III. 2 Primordial stars formed under the influence of a previous generation of stars, either by the ionizing or 
photodissociating radiation which they emit. 

Based on the enhanced cooling of the gas due to high H2 and HD fractions, it is expected that the typical mass 
scale of Pop III. 2 stars is significantly lower than that of Pop III.l stars (e.g. Uehara & Inutsuka 2000; Nakamura 
& Umemura 2002; Mackey et al. 2003; Machida et al. 2005; Nagakura & Omukai 2005; Johnson & Bromm 2006; 
Ripamonti 2007; Yoshida et al. 2007b). Following equation (4), the Jeans mass for gas that cools to the temperature of 
the CMB, which sets a rough upper limit for the mass of Pop III. 2 stars formed from partially ionized primordial gas, 
is 

where we have normalized to the same characteristic density at which the primordial gas fragments in Pop III. 1 star 
formation (e.g. Bromm & Larson 2004). At this fixed density, the Jeans mass is roughly an order of magnitude lower 
than expected for the case of Pop III.l star formation in minihalos. Also due to the lower temperature of the gas in the 
Pop III. 2 case, the rate of accretion onto a protostar is similarly lower (Yoshida et al. 2007b): 

3 

M^cc =^ 5 X lO-"^ (^^^ ' Mq yr-i , (22) 

compared to Mace ~ 10^^ M© yr^' for the case of Pop III. 1 stars (see equation 5). 

While enhanced molecule abundances are likely to result in lower characteristic stellar masses, other mitigating 
effects also come into play in the formation of second generation primordial stars. One factor which likely becomes 
important for shaping the stellar IMF, particularly in atomic cooling halos (e.g. Wise & Abel 2007b; Greif et al. 2008), 
is the development of supersonic turbulence (see e.g. Mac Low & Klessen 2004; Clark et al. 201 lb). Also, the degree 
to which the abundances of H2 and HD can be raised in the first galaxies is dependent on the strength of the molecule- 
dissociating radiation field generated by the first generations of stars (e.g. Wolcott-Green & Haiman 201 1). In the next 
Section, we shall see that an elevated radiation field may not only result in higher Pop III star masses, but may also 
result in the formation of the seeds of the first supermassive black holes. 



2.2 Suppression of Cooling by the Photodissociation of Molecules 

The assembly of the first galaxies becomes much more complex with the formation of the first stars, in part because 
they emit high energy radiation that alters the primordial gas in dramatic ways (e.g. Ciardi & Ferrara 2005). As in the 
case of primordial gas in the minihalos in which the first stars form, in the first galaxies one of the primary cooling 
processes is the emission of radiation from molecular hydrogen, and high energy radiation emitted by the first stars 



Before the adoption of this terminology, Pop III.2 was formerly referred to as Pop II.5 in the literature (e.g. Mackey et al. 2003; Johnson 
& Broinm 2006). 
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Fig. 6 The properties of simulated cosraologi- 
cal dark matter halos when the primordial gas 
first collapses to form Pop III stars, under the 
influence of vaiious levels of the background 
Lyman-Wemer photodissociating flux 
(shown here in units of erg s^' cm^^ Hz^' 
sr^'): star formation redshift Zcoll vs. /lw 
vs. Zcoii, Tvir vs. 7lw, and Mh vs. Jlw (clock- 
wise from top left panel). As the intensity of the 
background H2 -dissociating flux increases, a 
given halo must grow to a larger mass and have 
a higher virial temperature before the primor- 
dial gas can cool, collapse and form stars. At a 
flux of Jlw ^ 4 X 10"^^ erg s"' cm"^ Hz"' 
sr"' (i.e. 721 — 0.04), there is a steep increase 
in the rvi, and Mh of star-forming halos due 
to the rate of H2 photodissociation becoming 
comparable to the rate of H2 formation. From 
0'Shea& Norman (2008). 
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can easily destroy these molecules (e.g. Haiman et al. 1997; Omukai & Nishi 1999; Ciardi et al. 2000; Glover & Brand 
2001; Mackacek et al. 2001; Ricotti et al. 2001). So called Lyman-Werner (LW) photons, with energies 1 1.2 eV < hv 
< 13.6 eV excite H2, leading in turn to its dissociation into atomic hydroge ifl (Stecher & Williams 1967): 

H2-Fr^H5^2H. (23) 

With the destruction of H2 molecules, the primordial gas cools less rapidly and this signals a change in the rate 
at which gas can collapse into minihalos and form Pop 111 stars. An estimate of the minimum LW radiation field 
necessary to significantly delay star formation in a minihalo can be found by comparing the timescale ffomi for the 
formation of H2 to the timescale for its photodissociation. For a general radiation field the photodissociation time can 
be expressed as fdiss — 3 x 10^ ^21^ yr, where the specific intensity 7lw of the LW radiation field is defined as /lw 
= ^21 X 10" erg s"' cm" -2 Hz^ sr"' (e.g. Abel et al. 1997); here, ^21 

is a dimensionless parameter normalized to a 
typical level of the radiation field. To estimate the formation time we note that, as shown in the left panel of Figure 3, 
primordial gas collapsing into a minihalo is roughly adiabatic until its density rises to roughly « ~ 1 cm"^, at which 
point its temperature is T ~ 10^ K. Therefore, it is only at this characteristic density and higher that H2 is effective at 
cooling the gas, in turn leading to gravitational collapse and the formation of stars. It is the formation time of H2 in 
these conditions, which is ffonn — 10^ yr, that is to be compared to the photodissociation timescale fdiss- Equating these 
two timescales, we find a critical LW radiation field intensity of the order of ^21 — 10"^, at which the suppression of 
H2 formation and cooling slows the process of Pop 111 star formation in minihalos (see Kitayama et al. 2001; Yoshida 
et al. 2003; Mesinger et al. 2006; Wise & Abel 2007a; Johnson et al. 2008; Trenti & Stiavelli 2009). 

Figure 6 shows the results of cosmological simulations of the collapse of primordial gas into minihalos, under the 
influence of different levels of a constant LW radiation field. As the panels in the Figure show, for higher 7lw the 
primordial gas in a given minihalo collapses to form stars at lower redshift Zcoib when the halo has grown to a higher 
mass Myir and has a higher virial temperature Tyir. The results of these simulations corroborate our estimate of the 
critical LW background, as Tyir andMyij increase most dramatically at ^21 — 0.04. 



^ While for simplicity we limit our discussion to the photodissociation of Hi, HD molecules are also destroyed via this mechanism. 
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Fig. 7 The probability distribution of tiie 
Lyman- Werner flux 7lw to which atomic 
cooling halos ai'e exposed at ; = 10, for a 
model accounting for the clustering of the 
galaxies which emit LW radiation (black solid 
histogram) and for one in which galaxies 
are assumed to be distributed uniformly {red 
dotted histogram). The mean value of the LW 
background found in these models is /^i — 
40 (dashed blue line), which is well above the 
level needed to suppress the rate of Pop III star 
formation in minihalos (see Figure 6). While 
rare, a small fraction of halos are exposed to a 
LW flux J21 > 100, high enough to completely 
suppress star formation in atomic cooling halos 
until a > lO"* Mq black hole forms by direct 
collapse instead (see Shang et al. 2010). From 
Dijkstra et al. (2008). 
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Because the mean free path of LW photons is generally large, up to ~ 10 physical Mpc, a roughly uniform back- 
ground field is quickly established when the first stars begin emitting radiation (e.g. Haiman et al. 1997). We can 
estimate the level of the H2-dissociating background radiation, as a function of the cosmological average star forma- 
tion rate p* per unit comoving volume, by assuming that massive stars which live for a time produce the LW flux 
and that tJlw LW photons are produced for each baryon in stars (see Greif & Bromm 2006). We then obtain for the 
number density iiy of H2-dissociating photons 

nyC^riLW^ ^ 1+z , (24) 

'"H 

where mn is the mass of the hydrogen atom, ~ 0.76 is the fraction of baryonic mass in hydrogen, and the mass 
density in stars is ~ p*f*. Converting this to the photon energy density My = hviiy, we obtain an estimate of 7lw as a 
function of the star formation rate per comoving volume: 

Jlw - -r— = (1+^) ' (25) 

where c is the speed of light. In terms of 721, this is 

J.^o.2m]( .„_,..t:..^ ('-^y (26) 



V 10^ / VlO"^M0yr-iMpc-3 / V 10 

where we have assumed an average lifetime = 5 x 10^ yr for stars that produce the bulk of H2-dissociating radiation 
(e.g. Leitherer et al. 1999; Schaerer 2002). For a population of metal-enriched stars formed with a Salpeter-like IMF, 
as is inferred for the Milky Way today, T7lw — 4 x 10^; however, for metal-free stellar population with a top-heavy 
IMF, this can be as high as 77lw — 2 x 10^ (see e.g. Greif & Bromm 2006). In equation (26) we have normalized to 
an intermediate value, for simplicity; furthermore, while the star formation rate at very high redshift is not known, we 
have here normalized to a rough value expected in the standard ACDM picture of cosmological structure formation at 
z ~ 10 (see e.g. Tornatore et al. 2007; Haiman 2009; Trenti & StiavelU 2009). 

A further estimate of the cosmological background /lw, in particular that near the end of the epoch of reionization, 
can be found by assuming that the flux just above the Lyman limit (i.e. at hv > 13.6 eV) is sufficient to reionize the 
universe (see e.g. Bromm & Loeb 2003a; Shang et al. 2010) and that the sources producing the ionizing flux also 
produce a comparable flux in the LW energy range, 1 1.2 eV - 13.6 eV. Relating the number density of hydrogen nuclei 
at redshift z to the number density «y of ionizing photons required to keep hydrogen photoionized in the IGM, we 
obtain 
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R (pc) R (pc) 

Fig. 8 Properties of the primordial gas in a cosmological simulation, for various elevated levels of the LW background radiation field 
7lWj here assumed to be produced by stars with an effective surface temperature 10'* K, as is appropriate for metal-enriched Pop II stars. 
Clockwise from top left: number density, temperature, H2 fraction, and free electron fraction, each shown as a function of the distance from 
the center of an atomic cooling halo. A LW flux with 721 > 100 can prevent the formation of a sufficient fraction of H2 to cool the gas below 
~ 6 X 10^^ K, even at very high densities. Due to the large Jeans mass (equation 4) and accretion rate (equation 5) this implies, instead of 
hosting Pop III star formation, atomic cooling halos which are exposed to such an intense LW flux are likely to host the formation of a > 
lO'' M,.; black hole by direct collapse. From Shang et al. (2010). 



ny~Ny > (27) 

'"H 

where Ny is the number of ionizing photons per hydrogen nucleus required to keep the universe reionized and X2bPcrit 
is the cosmological average mass density of baryons at z = 0, expressed as a fraction £2\, of the critical density p^rit 
for a flat universe. Assuming that all LW photons which are emitted from sources within galaxies escape into the 
IGM, and taking it that only a fraction /esc of ionizing photons are able to escape due to the higher optical depth to 
photoionization, we find an estimate of the background flux as 

^^^^ 1 ^AVXWh(^^^^3^ (28) 

/esc 47r mu 

where again we have converted from photon energy density My = hvtiy to units of specific intensity as in equation (25). 
Expressing this in terms of 721, we have 



/.=^400(^ i^) , (29) 



12 



where we have normalized Ny to the value estimated by Wyithe & Loeb (2003), and /esc is normalized to a typical 
value found in cosmological radiative transfer simulations (e.g. Ricotti & Shull 2000; Ciardi & Ferrara 2005; Wise & 
Cen 2009; Razoumov & Sommer-Larsen 2010; Yajima et al. 201 1). 

This estimated level of the cosmological background radiation field during reionization is well above the critical 
level of 721 — 0.04 required for suppressing the rate of Pop 111 star formation in minihalos, and this may have important 
implications for the nature of the stars that are formed. In particular, under the influence of such an elevated LW 
background, due to the destruction of the H2 molecules which cool the gas, the temperature of the primordial gas 
when it finally collapses to form a star can be considerably higher than in the absence of a background H2-dissociating 
radiation field (O'Shea & Norman 2008). This, in turn, results in a higher Jeans mass and protostellar accretion rate, 
likely leading to more massive Pop 111 stars forming in the presence of a high LW background flux. 

While the LW radiation field is in general relatively uniform, near individual galaxies it can be locally higher than 
the cosmological average (see Dijkstra et al. 2008; Ahn et al. 2009), as shown in Figure 7. In rare regions where 
the LW background radiation is exceptionally high, a different outcome besides Pop 111 star formation in dark matter 
halos may result: the formation of a black hole by direct collapse (e.g. Bromm & Loeb 2003a). For this to occur, the 
LW radiation field must be at a level high enough to destroy molecules not just in the outskirsts of halos where the 
primordial gas begins to cool via emission from H2 molecules, but also high enough to destroy H2 even in the central 
dense regions of the halo (but see Begelman & Shlosman 2009; Mayer et al. 2010). 

Figure 8 shows the results of cosmological simulations from which the minimum /21 required for the formation of 
a black hole by direct collapse can be estimated. As shown in the bottom-right panel, for ^21 > 100 the H2 fraction 
in the gas is kept to a low level at which H? cooling does not lower the temperature of the gas significantly below 
the virial temperature of Tyir ~ 10^ K of the halo (Shang et al. 20100 Therefore, when the gas finally collapses, the 
accretion rate of primordial gas will be very high, of the order of ~ 0.1 M© yr^\ as can be seen from equation (5). 
This is roughly two orders of magnitude higher than the accretion rate onto Pop 111 protostars formed in H2-cooled gas 
at r ~ 200 K, and the result is predicted to be an extremely massive 'quasi-star' which quickly collapses to form a 
black hole with a mass > 10"* M© (e.g. Bromm & Loeb 2003a; Koushiappas et al. 2004; Begelman et al. 2006; Spaans 
& Silk 2006; Lodato & Natarajan 2006; Regan & Haehnelt 2009; et al. 2010b). While this level of the background LW 
radiation field is expected to be higher than the average, as shown in Figure 7, due to the clustering of the stars and 
galaxies producing LW radiation there may be a significant number density of black holes formed by direct collapse 
in the early universe. Indeed, some of these may be the seeds of the supermassive black holes observed at z < 6 (see 
e.g. Haiman 2009). 



2.3 The Impact of Radiation from Accreting Black Holes on the Primordial Gas 

In addition to the radiation emitted by the first generations of stars, black holes formed and assembled into the first 
galaxies can also produce radiation which dramatically impacts the primordial gas. In particular, the effects of the 
radiation emitted from black holes formed by direct collapse can be especially strong, as there is an initially large 
reservoir of gas that can be accreted onto the nascent black hole (see Johnson et al. 2011). To draw a comparison 
between the radiation emitted from stars in the first galaxies to that emitted during the accretion of gas onto black 
holes, we can calculate the temperature of the accretion disk and compare it to the typical effective temperature of 
a star. For a steady accretion flow, the temperature T of the accretion disk can be estimated by balancing the rate at 
which the disk is heated with the rate at which it cools. The heating is due to the gravitational potential energy of 
matter falling through the disk being dissipated by viscosity; for material falling through the disk and onto the black 
hole at a rate Mbh, the resultant heating rate per unit area F of the disk can be estimated on dimensional grounds as 
(e.g. Pringle 1981) 



It is important to note that the spectrum of the radiation producing the LW background must also be taken into account. While the results 
shown in Figure 8 are derived under the assumption that the LW background is generated by stars with an effective surface temperature of 
iff* K, appropriate for Pop II stars, higher levels of the LW flux are required to suppress H2 fonnation if, for instance, it is generated by 
massive Pop III stars with effective surface temperatures of ~ 10^ K (see e.g. Shang et al. 2010). 
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Assuming the disk is optically thick, then the rate at which the disk cools per unit area A can be estimated using the 
Stefan-Boltzmann law: A = GsqT'^, where cJsb the Stefan-Boltzmann constant. Equating these rates yields a temper- 
ature profile for the disk. The profile thus obtained is very close to the following formal solution, but for a correction 
near the inner edge of the disk r^n where the viscous heating rate goes to zero: 



T{r) 



3 GMbhMbh 



lO^K ( 



V lO^Mr: 



r 

10^ 



(31) 



In the second part of the equation we have normalized to a black hole mass of 10 M©, appropriate for the initial mass 
of a black hole formed by direct collapse. We have also normalized the radius to 10 Scharzschild radii = IGMqu/c^, 
which is well outside the inner edge of the accretion disk, < 3rs. Finally, we have assumed accretion to take place 
at the Eddington rate Mem, at which the outward force due to electron scattering of the emitted radiation balances the 
inward gravitational force acting on the accreting gas: 



47rGMBH'«H 



2 X lO"-"^ ( -- 
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/ Mbh 
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©yr 



(32) 



ecor vo.i/ VIO'^M©^ 

Here Or = 6.65 x 10^^^ cm' is the Thomson cross section for the scattering of photons off electrons, and e is the ratio 
of the radiated energy to the rest mass energy of the accreting material, normalized to a value appropriate for a slowly 
rotating black hole. 

As shown in equation (31), the temperature near the inner edge of the accretion disc of a rapidly accreting black 
hole can be as high as T ~ 10^ K. This is much higher than the effective temperature of even a very massive Pop 111 
star, which is roughly two orders of magnitude lower. In turn, this impUes that accreting black holes in the first galaxies 
emit both copious ionizing radiation and substantial LW radiation, as well as high energy X-rays (see e.g. Ricotti & 
Ostriker 2004; Kuhlen & Madau 2005). In the case of a black hole formed by direct collapse, the resultant photoheating 
of the gas in the host atomic cooling halo can drive its temperature to ~ 3 x 10"* K, as shown in Figure 9 (Johnson et 
al. 2011). Along with the associated high radiation pressure, this results in the expansion of the gas surrounding the 




Fig. 9 The projected number density {left panel), temperature {middle panel), and H II fraction {right panel) of the gas in the vicinity of an 
accreting 2.5 x lO'' M,.; black hole formed by direct collapse in an atomic cooling halo at z ~ 12. The ionizing radiation emitted from the 
accretion disk of the black hole has a strong impact on the gas in the halo, heating it to temperatures ~ 3 x 10'* K and causing expansion 
and disruption of the dense gas in the center of the halo from which the black hole feeds. From Johnson et al. (201 1). 
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black hole. The resultant drop in the density of the accreting gas translates into a decrease in the accretion rate of the 
black hole, which can be estimated by assuming gas which is gravitationally bound to the black hole falls towards it 
at the sound speed (see e.g. Bondi 1952). With the radius within which gas is bound to the black hole given by = 
IGMqh/ (cj + v|jj), where vbh is the velocity of the black hole relative to the gas, the accretion rate is estimated as 
the rate at which mass passes within a distance of the black hole: 

, , , , ' 47tG^Ml„umiin 

Mbh nrl^mun Kh + Cs) " = ^ , , , 

= 4 X 10- (^) (^^) f^Y' M, yr-' . (33) 

Vio^Mgy V0.6/ Vio2cm-3y Vio^Ky -"^ ^ ' 

In the second part of the equation we have assumed a black hole at rest with respect to the gas (vbh = 0) and we 
have again related the gas temperature to the sound speed using 'ikoTIl = pLmnc^ll. As the accretion rate is directly 
proportional to the density of the accreting gas and inversely proportional to its temperature, that the high energy 
radiation emitted from the accretion disk acts to heat and rarify the gas means that the accretion rate itself is regulated 
by the radiation generated in the process. Indeed, the Eddington rate given by equation (32) provides an estimate of the 
maximum rate at which gas can be accreted in the face of the intense radiation that is emitted. However, hydrodynamics 
calculations of accretion onto black holes formed in the first galaxies suggest that Mbh is on average well below the 
Eddington rate because of both strong radiative feedback during accretion (see e.g. Pelupessy et al. 2007; Alvarez et 
al. 2009; Milosavljevic et al. 2009; Park & Ricotti 2010; Johnson et al. 201 1) and low gas densities (e.g. Yoshida 2006; 
Johnson &. Bromm 2007). This poses a challenge for the rapid growth of black holes in the early universe. 

A further challenge to the model of black hole formation by direct collapse is the enrichment of the primordial gas 
with the first heavy elements (e.g. Omukai et al. 2008; Safranek-Shrader et al. 2010), which can easily cool the gas 
more efficiently than either H2 or HD molecules. We turn next to the broader question of how the first supemovae, 
which enrich the gas, transform the process of star formation in the first galaxies. 



3 Metal Enrichment and the Onset of Population II Star Formation 

We have seen that the characteristic mass of objects that form from the runaway gravitational collapse of gas, stars 
and in extreme cases black holes, depends critically on the temperature of the collapsing gas. The hotter the gas, the 
larger the Jeans mass and the higher the rate at which gas accretes onto the collapsed object. Therefore, the injection 
of heavy elements by the first supernovae represents a fundamental transition in star formation, in that new coolants 
are added to the primordial gas. As a result, the characteristic mass of stars formed from the first metal-enriched gas 
is likely to be lower than the characteristic mass of primordial stars. Here, we investigate the transition between these 
two modes of star formation. 



3.1 The First Supernovae and Metal Enrichment 

It is one of the hallmark predictions of modern cosmology that the first heavy elements, such as carbon, oxygen, and 
iron, are produced in the cores of stars and in supemovae, rather than in the Big Bang (e.g. Burbidge et al. 1957). 
Thus, when the first stars explode as supemovae the first metals, forged in their cores, are violently ejected into the 
primordial gas. In this, the first supemovae introduce not just new chemical elements, but also tremendous amounts of 
mechanical energy that dismpt their environments. Indeed, as given in Section 1, the definition that we have chosen 
for the first galaxies pertains to this: in the first galaxies, formed in haloes with virial temperatures Tvir ~ lO'* K, the 
gas can not be completely expelled by a single powerful supernova, as is the case in the minihalos hosting the first 
stars (see e.g. Bromm et al. 2003; Kitayama & Yoshida 2005; Greif et al. 2007; Whalen et al. 2008). 
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Fig. 10 The evolution of the remnant of a powerful Pop III.l supernova, exploding in a minihalo at z ~ 20. From left to right: the number 
density, temperature, and metallicity of the gas along the line of sight. From top to bottom: a time series showing the simulation 15, 100, 
and 300 Myr after the SN explosion. Each panel is 100 comoving kpc on a side, while the inlays show the central 10 comoving kpc. The 
metals are initially distributed by the bulk motion of the SN remnant, and later by turbulent motions induced by photoheating from other 
stars and the re-collapse of the shocked gas. As can be seen in the inlays, the gas collapses to high densities once again in an atomic cooling 
halo in which a first galaxy later forms from the enriched gas. From Greif et al. (2010). 

The effects of a powerful Pop 111.1 supernova on the primordial gas are shown in Figure 10, as gleaned from the 
cosmological simulation presented in Greif et al. (2010). Consistent with the results presented in Figure 1, the gas 
within the minihalo hosting the progenitor Pop 111 star is completely blown out into the surrounding IGM. There the 
primordial gas is shock-heated to several thousand Kelvin and enriched to metaUicities of up to ~ 10"^ within of 
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Fig. 11 The evolution of the remnant of a 
10^" erg Pop III. 1 supernova in a cosmological 
minihalo. The analytical solution {dashed 
lines) discussed in Section 3.1 accurately 
describes the expansion of the blast wave, the 
mass-weighted extent of which is shown by the 
black dots. As the results shown here are from 
a cosmological simulation, the inhomogeneous 
density field into which the shock propagates 
leads to some dispersion in the distance that it 
extends in different directions (gray region). As 
shown in Figure 10 for a similar cosmological 
simulation, after of the order of 10** yr the 
supernova shock stalls and the metal-enriched 
gas recollapses into the growing halo which 
hosted the progenitor star. From Greif et al. 
(2007). 
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the order of 10^ yr. The evolution of the supernova remnant can be well described analytically, as it passes through 
the four distinct phases of an explosion with energy /ssn = 10^^ erg in a medium with particle number density n < 1 
cm^^, as expected for a Pop 111.1 progenitor star with a mass of the order of 100 (e.g. Fryer et al. 2001; Heger 
& Woosley 2002; Whalen et al. 2008). At first, the blast wave from the supernova propagates outwards at a roughly 
constant velocity Vgh; in this, the so-called free expansion phase, the distance which the shock has traveled from 
the site of the explosion by time fjh is given simply by 

vshfsh ^ (^) %sh 3 " ' ( J^) PC . (34) 

At this stage all of the energy of the supernova is in the kinetic energy of the ejecta, which has an initial mass Mgj. 
When the shock has swept up an amount of mass comparable to the original ejecta mass, the shock enters the so-called 
Sedov-Taylor phase in which the energy of the blast wave is conserved while an increasing amount of mass Msw is 
swept up by the shock. In this phase we therefore have Vsh = '^'"sh/'^fsh — (2£'sn/^sw)'^^, which yields for the shock 
radius 

'"^'^ " ( ll%) ' ( l£^) ' ( ll^) ' ' ^''^ 

where we have used Msw = 47r/3rsh^/^mH«, with \i = 0.6, which is appropriate for an ionized primordial gas. The 
transition between the free expansion and Sedov-Taylor phase is evident in Figure 11, which charts the propagation 
of the blast wave of a powerful 10^^ erg primordial supernova in a simulated cosmological minihalo, similar to that 
shown in Figure 10. 

The third phase, also shown in Figure 11, sets in when a substantial fraction of the original energy in the blast 
wave has been radiated away, principally by recombination and resonance line cooling of the hydrogen and helium 
composing the primordial gas (Greif et al. 2007; Whalen et al. 2008), but also to some extent by bremsstrahlung and 
inverse Compton scattering of the CMB by free electrons, the latter being most important at high redshift due to the 
steep increase of the energy density of the CMB with redshift (e.g. Oh 2001). Known as the pressure-driven snowplow 
phase, at this stage the high pressure gas behind the blast wave powers its expansion, and the equation of motion thus 
becomes 
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Fig. 12 The metallicity of the gas enriched by 
the violent expulsion of metals in a powerful 
Pop III. 1 supernova, similar to that shown 
in Figures 10 and 11, as a function of the 
density and temperature of the gas. The highest 
metallicity gas is in the low density regions into 
which the blast wave propagates most rapidly. 
As also shown in Figure 10, the metallicity of 
the higher density gas re-collapsing into the 
halo hosting the progenitor star is Z ~ 10^"' 
Zq. From Wise & Abel (2008). 
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where F\, is the pressure in the hot bubble interior to the blast wave. As discussed in Greif et al. (2007), at the radius 
where the transition to the snowplow phase begins, the density profile of the gas is close to that of an isothermal 
gas, n 0= r^^; within this radius, the density profile is much flatter due to the strong photoheating of the gas by the 
progenitor star (e.g. Kitayama & Yoshida 2005; and Whalen et al. 2008). Therefore, as the pressurized bubble expands 

adiabatically, in the snowplow phase we have Msw °= ?"sh and F\, . This, in turn, allows a solution to the equation 

2/5 

of motion with r^h °^ f^]^ , just as in the previous Sedov-Taylor phase. 

The final transition occurs when the bubble behind the blast wave has cooled and the pressure behind the shock 
no longer affects it dynamically. While by this time a large fraction of the energy of the supernova has been radiated 
away, the momentum that has accumulated in the dense shell of gas that forms behind the shock is conserved. As the 



density profile of the ambient gas is still n < 



at this point, the conservation of momentum implies that the quantity 



Msv/Vsh °^ fsh dr^ih/dtsh is a constant. Thus, in this final phase of the supernova remnant rsh °^ tj^ , as shown in Figure 
11. 

The explosion of a Pop 111.1 star with a mass of ~ 200 M(;) as is shown Figures 10 and 1 1, is expected to release 
up to 10^^ erg as well as all of the up to ^ 100 in metals produced in the core of the star (e.g. Heger & Woosley 
2002; Heger et al. 2003; Karlsson et al. 201 1). The metal-enriched gas that is ejected into the IGM by the supernova 
explosion expands preferentially into low density regions, as shown in the middle row of panels in Figure 10. This 
can be seen more explicitly in Figure 12, which shows the metallicity distribution of the gas enriched by a similar 
primordial supernova as a function of density and temperature (Wise & Abel 2008). 

The dark matter in the halo hosting the progenitor star is not nearly so violently disrupted as is the gas swept up 
in the blast wave, and in fact the host halo continues growing until its gravity is strong enough for the cooling, metal- 
enriched gas to collapse into it again. As shown in the bottom panels in Figure 10, this occurs when the host halo has 
grown massive enough to host a first galaxy, as the gas is shock-heated to a temperature of 10"* K at the virial radius 
'"vir ^ 1 physical kpc from the center of the ^10^ Mq halo. As shown in both Figures 10 and 12, the metallicity of 
the gas that re-collapses into the growing host halo is typically ~ 10^^ Zm. Therefore, it is expected that stars formed 
in first galaxies enriched by powerful Pop 111.1 supemovae are likely enriched to this level (e.g. Karlsson et al. 2008; 
Greif et al. 2010; Wise et al. 2010). Such stars would be the first Pop 11 stars, and as we shall see many of these stars 
may still be present today, 13 Gyr after the formation of the first galaxies. 
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3.2 The Mixing of Metals with the Primordial Gas 



Here we consider two distinct situations in which the metal-enriched ejecta of primordial supemovae mix with the 
primordial gas, drawing on the results of the cosmological simluations of Pop 111 supernovae discussed in Section 3.1. 
Firstly, we shall estimate the timescale on which the primordial gas in the IGM that is swept up by the blast wave 
becomes mixed with the ejecta. Then we will turn to consider the likelihood that the primordial gas in minihalos that 
are overrun by the blast wave is mixed with the ejecta, thereby precluding Pop 111 star formation in those halos. 

When the supernova shock finally stalls after ^10^ yr, the dense shell of swept-up gas is accelerated towards 
the growing halo embedded in the underdense shocked gas. Such a configuration is Rayleigh-Taylor unstable and 
small perturbations of the shell can quickly grow, leading to mixing of the primordial gas in the shell with the metal- 
enriched gas in the interior As a stability analysis shows, a small perturbation on a length scale e << Zsh, where /sh is 
the thickness of the dense shell, will grow exponentially, at a rate 
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Here g is the acceleration of the dense shell in the direction of its interior, and and pb are the densities of the 
shell and the interior metal-enriched bubble, respectively. Assuming that Psh >> Pb, we can estimate the timescale on 
which the perturbation will grow as (e.g. Madau et al. 2001) 



dt 

where in the last expression the gravitational acceleration towards the growing host halo is taken to be g ~ GMh/r^jj, 
with Mh = 10^ M0, roughly the mass to which the host halo grows during the expansion of the blast wave. We have 
also used rsh = 3 kpc, which is roughly the spatial extent of the supernova shock when it finally stalls (e.g. Greif et 
al. 2007). Even for a shell as thick as ^ 100 pc, the timescale on which the metal-enriched interior material mixes 
with the 10^ M,7, of primordial gas swept up by the blast wave is much shorter than the timescale on which the gas 
re-collapses into the host halo, which is ~ 10^ yr Therefore, the gas which re-collapses into the host halo is expected 
to be well-mixed with the metals ejected in the supernova explosion. 

While it is thus apparent that the low density gas swept up in the IGM can be efficiently mixed with the metal- 
enriched material ejected in Pop 111 supernovae, the blast waves from these powerful explosions can also impact the 
more dense primordial gas inside neighboring minihalos. It is therefore another key question whether the metals are 
also mixed with this dense gas, as if so then when it collapses metal-enriched Pop 11 stars may form instead of Pop 111 
stars. In this case of a supernova blast wave overtaking a dense cloud of self-gravitating gas in a minihalo, there is the 
possibility of the dense gas becoming Kelvin-Helmholtz unstable, in which case vortices develop at the boundary with 
the fast-moving metal-enriched gas, and the two will mix with one another However, for this to occur the dense gas 
cloud must not be too tightly bound by gravity. In particular, for a given relative velocity between the minihalo and 
the blast wave, which we can take to be Vgh, the gas will be mixed due to the Kelvin-Helmholtz instability at the virial 
radius ryir of the halo, if (e.g. Murray et al. 1993; Gen & Riquelme 2008) 

where for the second expression we have used g = GM\ylr^^j^ and we have assumed a density contrast between the gas 
at the virial radius and that of the shock p^lpb = 10, consistent with the results of the cosmological simulations of 
Pop 111 supernvovae shown in Section 3.1. Thus, we see that the gas near the virial radius will be mixed if the speed 
of the shock exceeds the circular velocity Vcirc of the halo, which is likely the case for a minihalo with mass ~ 10^ 
M© at z < 20. 

However, as it is generally the dense gas embedded more deeply in the halo from which stars form, the metal- 
enriched gas may have to be propagating at a significantly higher velocity in order to impact the nature of star formation 
in the halo. To estimate how fast the shock must be in order to mix the gas a distance r from the center of the halo. 
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Fig. 13 The mixing of the primordial gas within a halo with mass Mh = lO' Mg, with metal-enriched gas overtaking the halo at a velocity 
Vsh = 30 km s^'. Shown is the metallicity of the gas at z = 9 (left panel) and later at z = 6 (right panel). The velocity field, along with 
contours of p = (solid lines), is depicted in each of the panels. The primordial gas is well-mixed with the metal-enriched gas at radii r 
> Jvi,-, and is mixed to a lesser degree at smaller radii, in basic agreement with equation (40). Note that the most dense gas in the center of 
the halo remains pristine; as it is from this gas that stars would most likely form, such a halo is likely to host Pop III star formation, despite 
being overtaken by high velocity metal-enriched gas. From Cen & Riquelme (2008). 

we can take it that the gas in the halo has a density profile that is roughly isothermal, with p ~ Pvir (/"//"vir)^- Using 
the same rough scaling also for the dark matter, we substitute p for Pvir and r for ryir in equation (39) to arrive at the 
following expression for the shock speed required for mixing via the Kelvin-Helmholtz instability: 

where we have implicitly assumed the same constant pb as in equation (39). Therefore, we see that it is only relatively 
fast shocks that are able to efficiently mix the metal-enriched material with the dense, pristine gas in the interior of a 
primordial minihalo. The results of this analysis are in basic agreement with the results of simulations of high velocity 
shocks impacting minihalos, as shown in Figure 13 (Cen & Riquelme 2008). 

To gauge the likelihood that a cosmological minihalo is indeed impacted by a Pop III supernova shock that is 
sufficiently strong to enrich the material in its central regions, we turn to Figure 14, which shows the average distance 
between Pop III star-forming minihalos in a cosmological simulation of the formation of a halo similar to that of 
the Milky Way (Gao et al. 2010). This Figure shows that due to the clustering of such halos, the average distance 
between them is smaller than would be expected from a simple estimate drawn from their abundance assuming a 
homogeneous distribution. In particular, especially at high redshifts (z > 20), the halos are closely clustered, with an 
average separation of roughly ^ 500 pc. From this we can estimate the average speed at which the blast wave from 
a Pop III supernova in one minihalo impacts its nearest neighboring minihalo, using the results presented in Figure 
11 . As can be seen from that Figure, the typical speeds with which the shock propagates at ^ 500 pc from the 
explosion site are roughly 20-40 km s^' . This is high enough to disrupt the gas near the virial radius of a neighboring 
halo, but not high enough to mix the metal-enriched ejecta with the dense star-forming gas in the interior of the halo at 
r < 0. 1 Tvir, as given by equation (40). Indeed, a similar result is found for the case of = 30 km s^ Mn simulations in 
which the mixing of the gas is resolved, as shown in Figure 13. Therefore, we conclude that the inefficiency of mixing 
poses a substantial challenge for the metals ejected in Pop III supernova explosions to enrich other star-forming halos 
and prevent Pop III star formation from occurring (see also Wyithe & Cen 2007; Wise & Abel 2008; Greif et al. 2010). 
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Fig. 14 Clustering of Pop III stai-forming minihalos in a simulation of the foimation of a Milky Way-like halo. Left panel: mean number 
of progenitor halo pairs, as a function of their physical separation at the three redshifts indicated in the panel. Right panel: mean separation 
of progenitor haloes in each of several similar simulations. Each halo has one neighbour within the distance given by the solid line, and 
10 neighbours within the distance given by the dashed lines. While the radiation produced by one Pop III star in a given halo is likely to 
impact its neighbors and perhaps delay star fonnation, as discussed in Section 2, it is in general unlikely that the metal-enriched ejecta from 
Pop III supemovae will mix with the dense star-forming gas in neighboring halos. From Gao et al. (2010). 



While here we have presented simple analytical estimates of the degree to which metals ejected in the first su- 
pernovae are mixed with the primordial gas via hydrodynamical instabilities, both in the IGM and in neighboring 
minihalos, other processes also contribute to mixing metals into the primordial gas (see e.g. Ferrara et al. 2000; Karls- 
son et al. 2011; Maio et al. 2011). Perhaps chief among these is the turbulence which develops as gas rapidly flows 
into the centers of the atomic cooling halos in which the first galaxies form (Wise & Abel 2007b; Greif et al. 2008) 
and acts to enhance the rate at which mixing takes place on small scales via diffusion (see e.g. Tenorio-Tagle 1996; 
Klessen & Lin 2003; Karlsson 2005; Pan & Scalo 2007). Once star formation begins in these halos, turbulent mixing 
is also facilitated by the energy injected by supernova explosions (e.g. Mori et al. 2002; Wada &. Venkatesan 2003; 
Vasiliev et al. 2008), and the fraction of un-enriched primordial gas in the first galaxies is expected to continually drop 
with time (e.g. de Avillez & Mac Low 2002). We turn next to discuss the impact that the first metals, once mixed into 
the primordial gas, have on the cooling of the gas and so on the nature of star formation. 



3.3 Metal Cooling in the First Galaxies 

In Section 2.1 we discussed how cooling by the molecule HD, which may be formed in abundance in partially ionized 
primordial gas, can lower the temperature of the gas to the lowest temperature possible via radiative cooling, that of the 
CMB. Here we draw on the same formalism introduced there to show how just a small amount of metals mixed into 
the primordial gas can allow it to cool to low temperatures even more efficiently. While a number of heavy elements 
contribute to the cooling of low-metallicity gas, here we shall take a simplified approach and focus only on cooling by 
carbon, which is likely to have been released in abundance in the first supernova explosions (e.g. Heger & Woosley 
2002, 2010; Tominaga et al. 2007). 

To begin, we note that once the first generations of stars form and a background radiation field is established, as 
discussed in Section 2.2, besides dissociating H2 this background radiation field can easily ionize neutral carbon (e.g. 
Bromm & Loeb 2003b). This makes available the potent coolant C II which, even in the Galaxy today, is important 
for cooling the gas to very low temperatures in dense star-forming clouds (e.g. Stabler & Palla 2004). To see how the 
presence of C II in the first galaxies affects the cooling of the gas, we first note that at low temperatures this ion is 
readily collisionally excited from its ground state to the first excited ^^3/2 state. The energy difference between 
these two fine-structure states is just 
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Fig. 15 Left panel: Cooling rates for primordial gas with a hydrogen number density iin = 1 cm^-* and the following fractions for the 
different species: Xhd = 10 , = 10^^ , Xjj+ = 10^'-' , X^- = lO^'* . The H2 cooling rate {long-dashed line) is plotted together with 
those of HD (solid), H-impact (dotted line) and e-impact Ho^ {short-dashed line). Right panel: Cooling rates of various metal species 
as a function of temperature, for a gas also with hydrogen number density «h = 1 cm^-*; for each metal species a fractional abundance of 
10^* is assumed. The cooling rate per C 11 ion is higher than the cooling rate per particle of any of the other species shown here, at Tg^s < 
100 K. From Maio et al. (2007). 



^^92K, (41) 

where here Vio denotes the frequency of the photon emitted in the radiative decay of the first excited state back to the 
ground state. As the energy difference is even smaller than that between the ground and first excited rotational states 
of HD, C II offers the potential to more efficiently cool the gas than HD, as even lower energy collisions are able to 
excite the ion. In addition, the Einstein coefficient for spontaneous radiative decay is Aio = 2.4 x 10^^ s^', almost 
two orders of magnitude higher than that for the 7=1—^0 transition of HD. 

We can obtain a conservative lower limit for the cooling rate of C II via this transition by considering the cooling of 
gas in the low density regime, in which the rate of collisional excitations is balanced by the rate of radiative decays, that 
is at densities n < «crit, where the critical density «crit is defined as that above which the rate of collisional deexcitations 
exceeds the rate of radiative deexcitations. For the transition of C II that we are considering, /icrft = 3 x 10^ cm^^^. In 
this case, the cooling rate is given as (e.g. Stahler & Palla 2004) 



Acii(« < «crit) - ^nonnYmhVioe^ = 1.5 x lO-^^ f (y^3^) ^ "^ergs-'cm-^ , (42) 

where Yw^h = 6 x 10^"* «h s^' is the rate at which a given C II ion in the ground state is excited due to a collision 
with a neutral hydrogen atom, gi is the statistical weight of the ;th excited state, Xqh is the fractional abundance of 
C II relative to hydrogen, and Tgas is the temperature of the gas. For densities « > iiait, in tuni^ the cooling rate varies 
linearly with the density of the gas, since in this case the level populations are in LTE, as given by equation (10), and 
the rate of radiative decay is no longer balanced by the rate of collisional excitation. In this case, we have 



Acii(« > »crit) - -«oAio/iVioe^ = 6 X 10'^^ 77^ ( " - J ^ '^^^ ergs"' cm-^ (43) 



The right panel of Figure 15 shows, along with the cooling rates of a number of other metal species, the cooling 
rate due to C 11 given above, as a function of gas temperature, for « < «crit- The cooling rates of oxygen, iron, and 
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Fig. 16 Temperature evolution of collapsing 
prestellar clouds with different metallicities, as 
a function of density. The evolutionary tracks 
of clouds with metallicities [Z/H]=-oo (Z=0), 
-5, -3, and -1 (-6, -4, -2, and 0) are shown by 
solid (dashed) lines. The dotted lines denote 
various values of the Jeans mass Mj. Near 
the critical density of H2, n ~ 10"* cm^^^, at 
which the cooling of primordial gas becomes 
less efficient, the temperature continues to 
drop as the gas collapses to higher densities 
for metallicities Z > 10^"* Z.,; this is due to 
the cooling provided by C II and other atomic 
species, as discussed in the text. At higher 
densities, n > lO'" cm^-*, the gas is able to 
cool efficiently even at a metallicity of Z ~ 
10^' Z.,, if dust is present, as is the case in the 
calculation shown here. From Omukai et al. 
(2005). 
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silicon that are shown can be obtained following the first part of equation (42) using the atomic data corresponding to 
those elements (see e.g. Santoro & Shull 2006; Maio et al. 2007). We note, however, that the cooling rate per C II ion 
is higher than that of any of the other metal species shown, as well as being at least an order of magnitude higher than 
the cooling rate per molecule of any of the primordial species shown in the left panel, at temperatures Tgas < 100 K. 
Therefore, we can focus on this chemical species as a means to derive a simple estimate of the minimum abundance 
of heavy elements required to significantly alter the cooling properties of the primordial gas, and perhaps thereby alter 
the nature of star formation. 

A rough estimate of the minimum carbon abundance required for the characteristic fragmentation mass to change 
from the relatively large value expected for primordial gas in the case of Pop III. 1 star formation can be found by 
considering the properties of the primordial gas when fragmentation takes place. At this stage, the so-called loitering 
phase in the collapse of the primordial gas in minihalos, Tgas ^ 200 K and n ^ 10'* cm^^ (e.g. Abel et al. 2002; Bromm 
et al. 2002). Hence, the Jeans mass (equation 4) is of the order of 100 M,7, and, if the gas does not cool efficiently then 
a massive Pop III star, or perhaps a binary or small multiple system, will likely form (e.g. Turk et al. 2009; Stacy et 
al. 2010; Clark et al. 2011a; Greif et al. 2011). However, if the gas cools to lower temperatures, then the Jeans mass 
becomes smaller and the gas is expected to fragment into smaller clumps; in turn, this is expected to translate into the 
formation of less massive stars. 

Following the discussion in Section 2. 1 , we note that in order for the gas to cool efficiently at this stage, the cooling 
rate must exceed the rate at which the gas is heated adiabatically by compression during its collapse (Bromm & Loeb 

2003b). Taking the adiabatic heating rate to be r^d ~ 1.5«A:B7"gas/fff, where % ~ (Gp)^^ is the free-fall time and the 
cooling rate Aqu is given by equation (43), this condition is satisfied if Xqh > 7 x 10^^. Assuming that all carbon is 
in the form of C II and taking it that the solar abundance of carbon is 3 x 10^"* by number, this yields a critical 
carbon abundance of [C/H]crit — -3.50 

While other elements, such as oxygen, iron, and silicon, also contribute to the cooling of metal-enriched gas, this 
abundance of carbon relative to the solar value is similar to what is found for the overall critical metallicity Zcnt / ~ 
10^^^ that is typically found in detailed calculations including atomic coolinJ^ (e.g. Bromm et al. 2001; Omukai et 
al. 2005; Smith & Sigurdsson 2007; Smith et al. 2009; Aykutalp & Spaans 201 1 ; but see Jappsen et al. 2009a,b). Figure 
16 shows the results of one such calculation, in which the temperature evolution of the gas is modeled as it collapses 
to high densities, for various values of the metallicity of the gas. For the case of metal-free gas, the temperature of the 
gas increases after the loitering phase at « ~ 10* cm^^; in this case, fragmentation at mass scales smaller than of the 



' Here we use the common notation for abundance ratios relative to those of the sun given by [X/Y] = logio(A'x/^Y) - logi()(^xWY)o, 
where A'x and Ny are the numbers of nuclei of elements X and Y, respectively. 

As the cooling rates of the various atomic species each contribute separately to the total cooling rate, it is the combination of their 
individual abundances which detennines whether the 'critical metallicity' is achieved (see e.g. Frebel et al. 2007). 
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Fig. 17 The state of the central gas cloud in a primordial minihalo at 100 yr (left panel), 300 yr (middle panel), and lO-* yr (right panel) 
after the formation of the first protostar in a very high resolution cosmological simulation. Shown here is the hydrogen density projected 
along the line of sight. Blackdots, crosses and stars denote protostars with masses below 1M|t,, between IM,;, and 3M,;,, and above SM r,, 
respectively. The gas fragments into a relatively rich cluster of protostars with a range of masses. Dynamical interactions can lead to the 
ejection of low-mass protostars, while more massive protostars tend to remain at the center of the cloud and continue to accrete from the 
suiTounding envelope of gas. Thus, while massive primordial stars are likely to form in such clusters, some fraction of Pop III stars with 
masses below ~ 1 M,.;, may also form. From Greif et al. (201 1). 

order of 100 M© is thus unlikely. However, when Z > lO^'* Z©, close the value we found above for the critical carbon 
abundance, the gas cools as it collapses to densities n > 10^ cm^^ and consequently the fragmentation scale decreases 
appreciably compared to the primordial case. Hence, less massive stars are likely to form in gas enriched to this level. 

We note also a second drop in the temperature of the gas at higher densities for even lower metallicities in Figure 
16; this decrease in temperature for Z > 10~^ Z© occurs because of dust cooling. While the dust fraction in extremely 
metal-poor gas is not known, dust formation in early supernovae (e.g. Nozawa et al. 2003; Schneider et al. 2004; 
Cherchneff &. Dwek 2010) may yield it high enough for the thermal evolution of the gas to be affected at n > 10^" 
cm^^, as shown here, even for such extremely low metallicities. In this case, the critical metallicity for low-mass star 
formation may be smaller than we estimated above for the case of cooling by atomic species such as C II, perhaps as 
low as Zcrit ~ IQ-^ Zq (see e.g. Schneider et al. 2006; Clark et al. 2008). 

In general, simulations of the evolution of low-metallicity star-forming gas give the same general result that the 
fragmentation scale, as well as the protostellar accretion rate, is higher for metal-free gas than for metal-enriched gas, 
and hence that the typical masses of Pop III stars are higher than those of Pop II stars (e.g. Bromm et al. 2001; Smith 
et al. 2007; but see Jappsen et al. 2009a,b). However, as shown in Figure 17, recent very high resolution cosmological 
simulations suggest that low-mass protostars formed in clusters may be ejected from the dense central regions of 
primordial minihalos due to dynamical interactions, in which case their growth may be limited due to the accretion of 
gas being dramatically slowed (Greif et al. 201 1; see also Clark et al. 201 la). In this event, it is possible that low-mass 
stars may indeed form from primordial gas, although they may only constitute a small fraction of all Pop III stars (e.g. 
Tumlinson et al. 2006; Madau et al. 2008). If their masses were less than ~ 0.8 M©, then such low-mass stars could be 
detectable as un-enriched dwarfs or red giants in the Galaxy even today (Johnson & Khochfar 2011), although there 
is a strong possibility that their surfaces would be enriched due to accretion of metals from the interstellar medium 
(see e.g. Suda et al. 2004; Frebel et al. 2009; Komiya et al. 2010). To date, however, no low-mass stars with overall 
metallicity below of the order of 10^^ Zq have been detected, which is consistent with the critical metallicity being 
set by cooling due to atomic species such as carbon and oxygen (see Frebel et al. 2007). 
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While the additional avenues for radiative cooling provided by even trace amounts of metals clearly alter the evo- 
lution of the gas and the process of star formation, other factors also play a role in dictating the thermal and dynamical 
state of the gas in the first galaxies. Magnetic fields may impede the large-scale collapse of the gas into dark matter 
halos (e.g. Schleicher et al. 2009; Rodrigues et al. 2010; de Souza et al. 2011) or alter the collapse of the gas at smaller 
scales during star formation (e.g. Kulsrud et al. 1997; Silk & Langer 2006; Xu et al. 2008; Schleicher et al. 2010a). 
Also, cosmic rays generated in the first supernova explosions are an additional source of ionization that can speed 
the formation of molecules and so enhance the cooling of the gas (see Vasiliev & Shchekinov 2006; Stacy & Bromm 
2007; Jasche et al. 2007). Finally, the impact of the turbulence generated by both the accretion of gas from the IGM and 
supernovae in the first galaxies may dramatically impact the process of star formation, in general acting to decrease 
the mass scale at which the gas fragments and forms stars (e.g. Padoan et al. 2007; Clark et al. 2008, 201 lb; Prieto et 
al. 2011). 



4 Observational Predictions and the Outlook for Identifying the First Galaxies 

While the enrichment of the primordial gas by metals ejected in the first supernovae is likely to preclude primordial star 
formation in a large fraction of the first galaxies (Johnson et al. 2008; Wise & Abel 2008; Greif et al. 2010; Maio et al. 
2010), it is also not likely that metal enrichment abruptly ends the epoch of Pop III star formation after the formation 
of the first stars. As discussed in Section 2.2, it is possible for the photodissociating background radiation established 
by early generations of stars to slow the collapse of the primordial gas, potentially delaying a large fraction of Pop III 
star formation and metal enrichment until later times. Also, as discussed in Section 3.2, the mixing of the first metals 
with the primordial gas, especially within minihalos, may not occur efficiently. Therefore, it is a distinct possibility 
that Pop III star formation continues well after the formation of the first stars (e.g. Scannapieco et al. 2003; Tornatore 
et al. 2007; Trenti et al. 2009; Maio et al. 2010), and that substantial primordial star formation may be detectable in 
the first galaxies. 



Fig. 18 The luminosity of Pop III star clusters 
in a first galaxy at j ~ 12, as a function 
of the time from their formation, in three 
recombination lines: Lycc (dot-dashed blue). 
Ha (solid red), and He II A 1640 (dashed 
black). The four panels correspond to four 
different combinations of stellar IMF and total 
stellar mass; these are, clockwise from top- 
left: twenty-five 100 M ., stars, two hundred 
fifty 100 M v, stars, one thousand 25 M ., 
stars, and one hundred 25 M,.) stars. The 
relative luminosities of detected He II and H I 
recombination lines can provide information 
about the stellar metallicity and IMF; more 
massive and more metal-poor stars emit more 
high energy radiation that can ionize He II, 
which leads to strong He II A 1640 emission 
relative to Ha and Lya. From Johnson et al. 
(2009). 
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Fig. 19 Stellar masses A^=^,niin of the lowest 
mass starburst observable through the detection 
of the Ha line (solid curves), the He II 11640 
line (dashed curves), or the continuum at 
1500 Mdash- dotted curves) with JWST, 
assuming an exposure time of 10^ s and a 
signal-to-noise ratio of S/N =10. Stellar 
masses derived from the Schaerer (2003) 
zero-metallicity starbursts with a standard 
Salpeter-like IMF, zero-metallicity starbursts 
with a top-heavy IMF, and low-metallicity 
starbursts, are shown in blue, red, and black, 
respectively. The right axis shows the masses 
^min = 10^ Mt.min of halos expected to host a 
starburst with stellai' mass min- The dotted 
curve shows the mass of a first galaxy-sized 
halo, with virial temperature Tjj^ = 10* K. 
From Pawlik et al. (2011). 
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It is therefore critical to predict observable signatures of Pop III star formation, in order that it can be identified in 
high redshift galaxies (e.g. Zackrisson et al. 201 1). Some distinctive signatures derive from the high surface tempera- 
tures of primordial stars, which arise due to a relatively low opacity in the stellar interior This low opacity translates 
into a smaller radii for primordial stars than for their metal-enriched counterparts. In turn, because stellar luminos- 
ity scales as °^ T^, for a given luminosity the surface temperature of a primordial star will be higher than a 
metal-enriched star (e.g. Siess et al. 2002; Lawlor et al. 2008). For very massive primordial stars, the surface temper- 
ature is very high, roughly ~ 10^ K (Bromm et al. 2001; Schaerer 2002). Owing to this high temperature, primordial 
stars emit copious high energy radiation, a relatively large fraction of which is able to ionize not only hydrogen (H l), 
but also helium (He I and He ll). 

As a substantial portion of the ionizing photons emitted from stars in early galaxies are absorbed by the relatively 
dense gas in the interstellar medium before escaping into the IGM (e.g. Wood & Loeb 2000; Gnedin et al. 2008; Wise 
& Cen 2009; Razoumov & Sommer-Larsen 2010; Paardekooper et al. 2011; Yajima et al. 2011), the energy in these 
photons is reprocessed into emission lines arising from the recombination of the ionized species (e.g. Osterbrock & 
Ferland 2006). For the case of primordial stars, because a relatively large fraction of the emitted radiation ionizes He II, 
the photons emitted during the recombination of He III to He II produce strong emission at characteristic wavelengths. 
The most prominent recombination line emitted from such He III regions, with a wavelength of 1640 A, emerges from 
the radiative decay of the lone electron in this ion from the n = 3 to the n = 2 stat^lH The most prominent emission 
lines from the recombination of ionized hydrogen in the H II regions surrounding primordial stars are the same as 
expected from metal-enriched stars, Lya and Ha, which arise from the radiative decay from the « = 2 — > 1 and « = 3 

2 energy levels of hydrogen, respectively. The key observational signature of primordial star formation, as opposed 
to metal-enriched star formation, is a relatively large ratio of the luminosity emitted in the helium line. He II A 1640, 
to that emitted in the hydrogen lines (see e.g. Tumlinson et al. 2001; Oh et al. 2001; Schaerer et al. 2003; Raiter et al. 
2010). 

Figure 18 shows the luminosity emitted in each of the three recombination lines mentioned above from an instan- 
taneous burst of Pop III star formation in a first galaxy formed in a halo of mass Mh 10^ at z ^ 12, as gleaned 
from cosmological radiative transfer simulations (Johnson et al. 2009). Each of the panels shows the Une luminosities 
for a different combination of the characteristic stellar mass of the stars (either 25 or 100 M©) and of the total stellar 
mass (either 2,500 or 25,000 M,:j). Even assuming such large characteristic masses for Pop III stars and that such a 
large fraction (either ~ 1 or ~ 10 percent) of the gas in the first galaxies is converted into stars, the luminosities of the 



' ' While photons are also emitted in transitions to the n = I state, the IGM is optically thick to these photons before reionization due to 
absorption by neutral hydrogen, and so they are not expected to be observable from the first galaxies. 
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recombination lines are likely to be too dim to detect with telescopes in the near future. To see this, we can estimate 
the total flux F that would be in these lines at z = as 



4lDj(z) 

where L is the luminosity in a given line and Dl is the luminosity distance to redshift z. At z > 10, even the most 
luminous line, Lya, would be seen at z = with a flux of of < 4 x 10^^" erg s^' cm^^, which is well below the flux 
limit of ~ 2 X 10^'^ erg s^' cm~^ of surveys planned for the JWST (Gardner et al. 2006; Windhorst et al. 2006). 

Instead of the first galaxies, hosted in halos with masses of ^ 10^ M,7, at z > 10, it thus appears likely that observa- 
tions in the next decade may reveal somewhat more developed galaxies hosted in more massive halos (e.g. Barkana & 
Loeb 2000; Ricotti et al. 2008; Johnson et al. 2009; Pawlik et al. 2011), although there is the possibility of detecting 
less developed galaxies if their flux is magnified by gravitational lensing (see e.g. Zackrisson 2011). As shown in 
Figure 19, the JWST is predicted to be capable of detecting both He II A 1640 and Ha from metal-free starbursts in 
halos with masses > 3 x 10** M©, if the IMF is very top-heavy. Though, if the typical mass of Pop 111 stars is < 50 M©, 
it is likely that He II A 1640 will only be detectable from significantly more massive stellar clusters, expected to form 
in similarly more massive halos. However, because more massive halos are formed from the mergers of smaller halos 
which themselves may have hosted star formation, it may be predominantly metal-enriched Pop 11 stars that form in 
the galaxies which will be detected by the JWST (e.g. Johnson et al. 2008)M 

That said, there is the possibility that substantial Pop 111 star formation takes place even well after the epoch of 
the first galaxies (i.e. at z < 10), either due to inefficient mixing of primordial and metal -enriched gas (e.g. Jimenez 
& Haiman 2006; Pan & Scalo 2007; Wyithe & Cen 2007; Dijkstra & Wyithe 2007; Cen 2010) or to the collapse of 
primordial gas into late-forming atomic cooling halos (e.g. Tornatore et al. 2007; Trenti et al. 2009; Johnson 2010). 
Pop 111 star formation at such late times could be detected more easily, in large part because the emission line flux 
increases strongly with decreasing redshift, as shown in equation (44). However, at such low redshifts the background 
ionizing radiation field that builds up during reionization can strongly inhibit the infall of primordial gas into halos, 
limiting the amount of Pop 111 star formation that can occur even in metal-free galaxies (e.g. Efstathiou 1992; Gnedin 
2000; Tassis et al. 2003; Dijkstra et al. 2004). 



5 Summary and Conclusion 

In this Chapter we have discussed a wide range of the physical processes that must be accounted for in the theoretical 
modeling of the first galaxies. We have made important distinctions between the formation of the first stars in minihalos 
and star formation in the atomic cooling halos hosting the first galaxies, highlighting how the cooling properties of the 
gas assembled into the first galaxies are altered by high energy radiation and by the injection of heavy elements from 
the first supemovae. While this has not been a complete review of the theory of the formation of the first galaxies, it 
has hopefully served to illustrate, from basic principles where possible, much of the physics that comes into play in 
their study. The reader is referred to the many excellent articles in the bibliography below for more in-depth study on 
the topic. 

In closing, it is critical to point out that without accounting for all of the effects we have discussed together, one is 
left with an incomplete understanding of the first galaxies. For instance, as we have seen, the radiation from the first 
stars can ionize the gas and trigger HD cooling, but it can also easily destroy H2 and HD molecules. As well, while the 
first supemovae may enrich much of the gas from which the first galaxies form to a level above the critical metallicity 
needed for low-mass Pop II star formation, much of the dense gas in minihalos may not be efficiently mixed with 



It is also likely that other, complementary next generation facilities, such as the Atacama Large Millimeter Array (e.g. Combes 2010), 
will detect only metal-enriched star-forming galaxies. 
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the metal-enriched ejecta. Similariy, black holes may only form by direct collapse in rare regions in which the LW 
background radiation field is elevated, but the same stars which likely produce this radiation may also enrich the gas 
when they explode as supernovae, possibly precluding this mode of black hole formation. A complete and consistent 
picture of the formation of the first galaxies only emerges when accounting for star and black hole formation, metal 
enrichment, and radiative feedback all together in the full cosmological context. 

Making this task especially daunting is the range of scales that must be taken into account. The gas clouds which 
collapse to form stars are on sub-parsec scales, metal-enrichment from the first supernovae occurs on parsec to kilo- 
parsec scales, and the radiation emitted by the first stars can impact regions on kiloparsec or even megaparsec scales. 
Thus, simulations must ultimately resolve an enormous range of scales in order to capture all of the important physical 
processes that come into play. While we have introduced the results of numerous analytical calculations and simula- 
tions, none of them alone captures all of the processes we have discussed simultaneously. Indeed, this stands as one of 
the primary challenges to making detailed predictions of the nature of the first galaxies. 
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